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We use quasiparticle description of deconfined matter in nuclear collisions at finite temperature 
and chemical potential. We assume that evolution of expanding system is isentropic. Using theoret- 
ical formulas and experimental meaning for the average multiplicity of charged and neutral particles 
and also of net nucleons (TV — N) per unit central rapidity, we calculate the initial temperature 
To, and the volume Vb for the collisions of heavy nuclei (Pb+Pb and Au+Au) at SPS and RHIC 
energies. At calculations we use the conservation of entropy and of number net nucleons in initial 
plasma stage and on the stage of phase transition at temperature T c , where we take into account 
the phase of constituent quarks m q and m s . We use at calculation the effective quasiparticle model 
with decrease of thermal gluon mass m g (T) at T — > T c from above. This model agrees with lattice 
data. The particle ratios for SPS and RHIC are defined by early chemical freeze-out at tempera- 
ture close to T c (at values of chemical potential [is = 247 MeV and 50 MeV for SPS and RHIC (at 
^J~s~ — 130 GeV) correspondingly). The spectra of baryons and mesons are calculated at temperature 
of thermal freeze-out Tf — 120 MeV. The quantitative characteristic of these spectra (for example, 
the normalization for baryons and mesons) are determined in the present model by initial state - 
quark-gluon plasma. At calculation of spectra we assume that averaged transverse flow velocity v± 
have both direct particles and paternal resonance. We show that values of v± increase for energy 
RHIC in comparison with SPS. We had investigated also the analogous problem of nuclear collisions 
in quasiparticle model with phenomenological parametrization of running coupling G(T), which 
increase at T — » T c from above (in difference from effective quasiparticle model). We have here 
agreement with new lattice data. We do not found in this model quantitative difference for spectra 
in comparison with effective quasiparticle model. However it can be shown, that in this model we 
have dramatic situation for energy loss of gluon jets in plasma. 

PACS numbers: 12.38.Mh, 24.85.+p, 25.75.-q 



I. INTRODUCTION 

There is considerable interest to heavy ion reaction 
at sufficiently high energies, because they provide the 
possibility for conditions of the transitions to deconfined 
state of hadronic matter-quark-gluon plasma (QGP). At 
Pb+Pb and Au+Au collisions there is considerable num- 
ber of net nucleons (N — N) in central interval of rapidity, 
which it is necessary take into account on plasma and 
hadron stage. Therefore we have the baryon chemical 
potential \lb + 1 0. 

In the works a number of authors (P.Braun-Munzinger 
et.al. 0) was shown, that statistical thermal model gives 
good description of particle ratios in central collisions 
of heavy nuclei at high energies (for SPS and RHIC). 
The temperature T for all this is close to critical tem- 
perature T c ps 170 MeV obtained from lattice Monte- 
Carlo simulations of QCD, and fiB is close to 50 MeV 
for RHIC(at -^5=130 GeV) This correspond to chemical 
freeze-out, when inelastic collisions cease to be impor- 
tant and the particle composition is fixed. The trajec- 
tory of chemical freeze-out T(/j,b) was obtained also in 
the work Q on the basic of relation for hadronic gas 
< E h > I < N h >= const ps 1 GeV. For SPS was found 
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the value fis ~ 250 — 260 MeV. One can to suppose that 
trajectory T(/i#) for various energies of nuclei is close to 
boundary of deconfined quark-gluon plasma phase. For 
all this the potential /i s correspond to strangeness neu- 
trality (i.e. overall strangeness = 0). This correspond 
in fact to hadronic part of mixed phase, where there are 
strange mesons and baryons (i.e. /i s ^ 0). 

In the present work we try to define the physical char- 
acteristic of initial plasma phase To and Vq. For this aim 
we use hypothesis of conservation of entropy and of num- 
ber of net nucleons in initial stage (where fi s = 0) and on 
stage of phase transition at T = T c , where we take into 
account also of constituent quarks phase. 

The plasma at finite temperature was considered per- 
turbatively up to order 0(g 5 ) Q. However in the ex- 
perimentally accessible region (close to temperature T c 
of phase transition) the strong coupling constant is suf- 
ficiently large: g ~ 2. The perturbative expansion in 
powers of g probably gives bad convergence. Recently in 
series of papers 0, la 01 f° r description of quark-gluon 
plasma have been used quasiparticle model, which allow 
apparently to sum up partially the perturbative expan- 
sion. At high temperature the plasma consist of quasifree 
quarks and gluons and probably is described well by per- 
turbative theory. The lattice model show rapid increase 
of entropy density and other values near by T c in SU(3) 
gluon plasma. In quasiparticle model the interacting 
plasma of quarks and gluons is described as system of 
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massive quasiparticles. 

For thermal momenta k ~ T the quark-particle exci- 
tation and transversal gluons will propagate with disper- 
sion relations 

uf(k) * rn 2 (T) + k 2 , (1) 
m?(T) = mli + Ih, 

where 11^ are given by the asymptotic values of the hard 
thermal self energies 

N 2 — 1 u 2 

n g = a^mo + w,), ^ 2 = i^[ r2 + 3 G2 ' ( 2 ) 

The contributions 11^ are generated dynamically by the 
interaction within the medium, and TOj(T) acts as effec- 
tive mass. The quantity G 2 (T) is to be considered as 
effective coupling in quasiparticles model. The entropy 
and energy density for example for gluon plasma takes 
the form: 

oo 

do f o 4k 2 + 3m 2 (T) 



oo 

e g (T) = ^ J k 2 f{E k ){E k )dk + B(T), (4) 
o 

where f(E k ) = [exp(f^) - E k = ^Jk 2 +m 2 g (T) 

The value B(T), introduced in expression for energy 
(and also for pressure) is necessary in order to maintain 
thermodynamic consistency, that is must be fulfilled the 
relations e + p = sT and s — The function B(T) 
compensates the additional T derivatives from m(T) in 
pressure and it represent in integral form in the paper Q . 
In our paper we do not use the function B(T). 

However the phenomenological parameterization of 
coupling constant G 2 in accordance with perturbative 
QCD in quasiparticle model 0, EJ is very likely un- 
just close to temperature of phase transition T c , where 
it should be expected of nonperturbative dynamics. In 
this model we have the increase of thermal gluon mass in 
vicinity of phase transition. 

In Sec. [H] we give the brief review of phenomenologi- 
cal model of confinement ^(Ji which agrees with SU(3) 
lattice data. This model is applied also to system with 
dynamical quarks. 

In Sec. lIIII we calculate for isentropic evolution of quasi- 
particles system with chemical potential /is ^ the var- 
ious physical characteristic: the initial temperature To, 
the volume Vo, the initial entropy So etc. for nuclear 
collisions at SPS and RIIIC energies. We use the values 
T c and [Ib close to above-mentioned: T c = 170 MeV and 
He = 50 MeV and 247 MeV for RHIC and SPS. At cal- 
culation we use also additional information - experimen- 
tal meaning of average number of net nucleons in central 



region of rapidity for SPS and RHIC. In result we have 
obtained the equations (|2*2l , l(2lfl) for definition of T . Si- 
multaneously must be fulfilled the relations (|24|l .which 
follow from conservation of entropy and of number net 
nucleons. 

Besides of study of characteristic of initial phase we de- 
scribe in Sec. lIIII another new aspect - we show, that from 
conservation of entropy and of number of net nucleons 
follow, that massive constituent quarks [m q and m s ) ap- 
pears with decrease of number of degrees of freedom (i.e. 
with effective number) in the presence of octet of pseu- 
dogoldstone states. We show also, that with the same 
effective number of degrees of freedom appears hadrons 
and resonances in hadron part of mixed phase (for SPS 
and RHIC). 

Certainly it should be noted, that in nucleus- nucleus 
collisions at high energies there is possibility of two stages 
of equilibrium. The gluons reach of thermal equilibrium 
to a considerable extent faster, than quarks. According 
to estimations, equilibrium time of gluons is r g ~ 0.3 
fm/c and for quarks one is r q ~ 2 fm/c, i.e. it is possible 
the production of hot glue at first stage • The corre- 
sponding estimations of initial conditions for hot glue in 
effective quasiparticle model there are in the work 

In Sec. lIIII we investigate the initial condition (To, Vo) 
and evolution of second plasma stage, where is the equi- 
librium also for quarks. In this paper we do not consider 
the particle ratios. We have considered the ratios for var- 
ious baryons and mesons with accounting of resonance 
decays in the work [33l |. The particle ratios in the main 
agrees well with experimental data for SPS and RHIC 
(Vs = 130 MeV. 

In Sec. lIVI we calculate in effective quasiparticle model 
the some baryons and mesons (pions and kaons) spec- 
tra with accounting of resonance decays and transverse 
flow. We use for heavy nuclei the thermal freeze-out tem- 
perature Tf — 120 MeV. The normalization of baryons 
and mesons spectra here is no free parameter, but is de- 
fined by initial condition in plasma - by initial entropy So 
and by condition of entropy conservation S/ = Sq. We 
show that values of averaged transverse flow velocity v± 
increase for RHIC in comparison with SPS. 

The main object of present work is no study the whole 
of spectra for SPS and RHIC, but study of dependence 
spectra from initial state and from stage of phase transi- 
tion plasma into hadrons. 

Therefore in Sec. [V] we consider the analogous isen- 
tropic problem of nuclear collisions by the use of ordinary 
perturbative theory up to order 0(a s (T)) with QCD pa- 
rameter A ~ 0.2 GeV. We calculate here also the values 
T , Vo, So- We do not find noticeable difference for spec- 
tra of particles in comparison with effective quasiparticle 
model. However this model disagrees with SU(3) lattice 
data in region of phase transition (i.e. near of temper- 
ature T c ). In Sec. we consider also the quasiparticle 
model with phenomenological paramctrization of running 
coupling G(T), which gives good fit of new lattice data. 
We calculate here also the initial parameteres Tq, Vq,So- 
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The coupling strength Gt here increase at T — > T c from 
above (unlike effective quasiparticle model). However 
baryons and mesons spectra in this model practically do 
not differ from spectra in effective quasiparticle model. 
Thus the spectra weakly depend on region of phase tran- 
sition. 

However it can be shown that in ordinary petur- 
bative model (taking into account the hot glue) we have 
too great energy loss AE g of gluon jet at RHIC energy: 
~ 80 — 90% (unlike effective quasiparticle model) . Thus 
the jet quenching essentially depend on description of 
phase transition stage. We show also, that in second 
above-mentioned model the situation with energy loss 
AE g is more dramatic. 

In Sec. lVIl - conclusion. 



II. EFFECTIVE QUASIPARTICLE MODEL 

In the papers 0, for gluon plasma was used the 
phenomenological parameterization of coupling constant 
G 2 (T) in accordance with perturbative QCD with the 
two fit parameters T s and A. The resulting equation of 
state (EOS) is in good agreement with lattice data over 
a wide temperature range between T c and 5T C . However 
in spite of this the thermal gluon mass have strong in- 
crease in the vicinity of the phase transition at decrease of 
T. Furthermore, close to a phase transition the coupling 
constant G S {T C ) ~ 2, that is correctness of perturbative 
of calculations is questionable. 

Further it is expected that correlation length £(T), 
which is proportional m^ 1 , grows when T — ► T c from 
above (mo is gluonic Debye mass). The lattice calcula- 
tions shows that mu drops by factor of 10 when T de- 
crease from 2 T c to T c p| Also for three colors the value 
£(T C ) remains large, but finite. In hard thermal loop 
(HTL) perturbation theory the thermal mass m g (Tland 
gluonic Debye mass mn are connected by relation 9] 



m D 



(5) 



However in nonperturbative theory this relation can be 
broken down 

Recently in the paper [Tcj was considered interesting 
phenomenological model of confinement. In region below 
T c in a pure SU (3) gauge theory are color singlet of heavy 
glueballs. Approaching T c the gluons are liberated, fol- 
lowed by sudden increase of energy and entropy density. 
When approaching the phase transition from above, the 
decrease of thermodynamic values is not caused by in- 
crease of masses (unlike 0, @), but caused by decrease 
of number active degrees of freedom. When T comes 
closer to T c , more and more of gluons are absorbed by 
glueballs. It is assumed that thermal gluons mass m g (T) 
follows roughly the behavior of the Debye mass, i.e. it 
decreases. However the entropy density s(T) will over- 
shoot the lattice entropy because light masses near T c . 
This difference may be accounted in quasiparticle model 



by modifying of number effective degrees of freedom in 



g g -> C(T)g g 



(G) 



The explicit value C(T) may be estimated as the ratio 
of the lattice entropy and entropy © with a dropping 
mass m g (T). At T ^> T c we expect C(T) ~ I and near 
T ~ T c we have G(T) < 1. At T < T c it can esti- 
mate C(T C ) ~ 0.2 from lattice data. The value C(T) is a 
smooth, monotonously increasing function with T. 

It is assumed also that thermal mass (for example glu- 
ons mass) have form m g (T) = G(T)T. The lattice results 
for Debye mass mjj (T) can be parameterized well by for- 
mula 



1/2 



m D {T) ~ Gl T(l + *-^= G T{l + 5-^f 

(7) 

where Gi~1.3,/3~0.1 and there is small gap at T = T c 
with 5 ~ 10~ 6 . The value G\ is determined here by the 
asymptotic value of thermal gluon mass, which coincide 
with lattice mass for instance at T = 3T C . Below of 
this value the explicit HTL resummation is expected to 
fail [l3|. One can to use the asymptotic value G(T) with 
renormalization group inspired parametrization 



G 2 (T ;A1 = 0) = 



24tt 2 



(lliV c - 2N f ) ln(A(T + T s )/T c ) 



(8) 



The parameteres T s and A are used for fitting new lat- 
tice results with N f =2 and N/=4 at T = 3T C (and 
also for Nf=3 6J). The thermal gluon mass at chemi- 
cal potential n q —0 is analogous with formula (JTJl where 
Gi -> y/N c /& + N f /12G . Effective coupling G have 
the form Gn = , 90 The gluon and quark mass 

y/llN c -2N f ^ 1 

coincide with lattice data at go ~ 9.4 for Nf—2 and at 
go ^ 9.8 for N f =3. That gives G ~ 1.886 for N f =3. 
The like value Go =1.9 we have in formula (7J. The such 
value Go — 1.9 we use in calculations. It should be noted, 
that lattice results at finite chemical potential [i q point to 
greatly weak dependence of asymptotic thermal masses 
from /j, q [l4l| . 

Since (in accord with assumption) the values m g (T) 
and tod(T) have similar trends near T ~ T c , the mass 
m g (T) is parameterized by analogy with formula Q. In 
order to account for certainties, were investigated 01 
some a range of values Go, /3 and S. 

A decreasing effective coupling strength G(T) at 
T — > T c from above can be understood, since at decreas- 
ing of T more and more gluons become confined and form 
heavy glueballs. The effective glueball exchange interac- 
tion between gluons are reduced. The total interaction 
can be interpreted as a superposition nonperturbative 
multigluon and weak glueball exchange. The good fit [T^j 
for C{T) (JSJ) have form similar with (J7J) 



C(T,T c ) = C Q [l + 6 c -j, 



/3c 



(9) 
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where C ~ 1.25, S c ~ 0.0026, & ~ 0.31. 

The relations J7J), © give good description of SU(3) 
lattice data for e/T 4 , 3s/4T 3 , 3p/T 4 for T close to T c . 

Further it is possible to extrapolate effective quasi- 
particle model to system with dynamical quarks. Un- 
fortunately no lattice data on thermal masses with dy- 
namical quarks are available. We assume that N c and 
Nf dependence of m g and m q are given by formu- 
las ©, where effective coupling G(T) have the form: 
G(T) = G (l + 5- T c /T f. For example, the thermal 
masses of u and d quarks are: 



where 



N 2 



+ 2m q0 G(T) ] 



16N C 



-(T 2 



N*-l 
)G 2 {T). 



{T 2 + 4) + (10) 



For s quarks we replace m q o — > m s o, fi q — > 0. It is 
assumed that parameterization of function C(T, T c ) in 
the presence of quarks and parameterization for glu- 
ons are similar, with some variation of parameters. For 
examplejp = 1-25, S c ~ 0.02, (3 C ~ 0.28 for two light 
quarks (lp|. The relations of type 10 for C(T,T C ) and 
relations of type (|10|) for m q and m g are inserted into for- 
mulas for entropy density of quasiparticles. For q(u, d) 
quarks we have the entropy density: 



s q {G Q ,m q ,fi q ,T,T c ) = 
= s ql (G ,m q ,fj, q ,T,T c ) + 
+ s q i(Go,m q ,-n q ,T, T c ) - 

_ t2_9q (Go,m q , fi q ,T, T c ) 

- Fi(G ,m q ,-fiq,T,T c )) 



(11) 



where 



Sql(G ,m q , fl q ,T,T c ) = 



_T 3 9« 



2tt 2 



dx- 



(4x 4 + Sx 2 ^ 1 )C(T,T c ) 



3^ + 4^[eV^^-^ + l] 



The formula l|ll|) gives correct term already in order g 2 
at perturbative expansion. 

For s quarks m q — > m s , and /i s = 0. For gluons the 
entropy density is: 



s g (G ,m g ,T,T c ) 

T 3 g g 



(12) 



2tt 2 



J ^ (4x 4 + 3. 2 ^P)C(T,T c ) 
o zJx 2 + ^±{J x2+ ^ 1) 



For the net nucleons density we have: 



ni(G ,m q ,fj,q,T,T c ) = (13) 
-jrp9q[Fi{Go,m q , fj, q ,T,T c ) - Fi(Gq, m q , —fj, q ,T, T c )], 



Fi(Go,m q ,^ q ,T,T c ) = j dx- 



x 2 C{T,T c ) 



We shall show below that in massive constituent quarks 
phase and hadron part of mixed phase we have also de- 
crease of the number of effective degrees of freedom. 



III. INVESTIGATION OF PHYSICAL 
CHARACTERISTIC OF INITIAL AND MIXED 
PHASES 

In "standard model" of relativistic nuclear collisions 
after a pre-equilibrium period the system becomes ther- 
malised in form of quark-gluon plasma with formation 
time To- The system then expand, cools off and at the 
time t c it begins to hadronize. At temperature T c ap- 
pears the mixed phase — at first plasma part at r c , then 
hadronic part of mixed phase at tjj- At temperature 
near to T c ~ 170 MeV particles first undergo a chemical 
freeze out, and then a thermal freeze out at temperature 
T f ~ 100 - 140 MeV. The particles ratio is defined by 
temperature close to T c . For A — A collisions we have 
the initial volume Vq = ttR 2 a t , where Ra — 1.2 A 1 / 3 . 
We consider isentropic evolution. The entropy per unit 
rapidity is conserved during the expansion [lg, H3 • The 

dS 
dy 



value ^| = ttR\s(t)t — const(^-) y=0 determine ini- 
tial entropy density, where ^ is the average multiplicity 
charged plus neutral particles per unit central rapidity. 
In this case part of initial energy goes into collective mo- 
tion of the expanding system, that is the initial energy 
density must have been higher than in the free flow sce- 
nario. In the work was investigated the isentropic 
longitudinal expansion of ideal gas of gluons and three 
flavors of massless quarks. In this work was found the 
dependence eo — (^) 4 ^ 3 f° r initial energy density in 
plasma. In this paper we have investigated more common 
case for quasiparticles with fig ^ 0. The quasiparticles 
density in plasma phase is: 

(14) 



n 



where n qi n Sl n g is density of u,d,s quarks and gluons. 
We have here: 



n q (G a ,m q , ^ 9 ,T T C ) 



(15) 



2tt 2 



[Fi(G , m q , T , T c ) + Fi(G , m q , -/j, q , T , T c 



n s (Go,m s ,T ,T c ) = / dx — , == , (16) 



n g (Go, m g ,T , T c ) 



2tt 2 



dx- 



+ 1 



x 2 C(T Q ,T c ) 



(17) 
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By analogy it can be written for summary entropy den- 
sity: 



so 



(18) 



where s q , s s ,s g given by (fTTJl . fT^I . We must express the 
initial temperature To across the number of secondary 
particles per unit central rapidity. For convenience of 
calculations we introduce the designation: no = TqHqi, 
,s = TqSqi, £o = Tq€oi, that is we single out the fac- 
tors Tq, Tq. We single out also the constant factors 
uqi = 12n 02 /(27r 2 ) , e i = 12e) 2 / (27r 2 ) Thus we have 
To = (no/noi) 1 / 3 , e = (n /n 01 ) 4 / 3 12e 02 /(27r 2 ). Hence 

we have here e = (jfc) /, where / = t02 f^J)}f 3 ^ ■ 

Thus the initial temperature To have been expressed 
across the number of secondary per unit central rapidity: 



To = 



/ dN 



\dyV 



(— Y 



(19) 



For A 
the form p]j 



A collisions the value 4rr can be extrapolated by 



AA 



dy 



d V Jpp 



(20) 



For pp collisions we have good approximation per unit 
central rapidity pcj : 



dN 
dy 



— ~0.81n- 



y=o 



(21) 



The value a describes the amount of rescattering in 
the interaction of the nuclei. From p — A data it can 
expect a ~ 1.1. The value a = 1.1 describes the result 
for various light ion quite well 0|- The data for SPS 
(V« = 17.2 GeV) is a = 1.08 ±0.06 (i.e. close to (p - 
A) data). We have ^ ~ 803 (for a =1.1) The data 

of Phenix collaboration show increase of ^ from SPS 
(Na49) to RHIC (y/s = 130 GeV) on value ~ 1.7, i.e. 
on value ^"(Vf-^o) p or gpg we assum e a = 1.1. From 

ln(ys=17.2) 

formula J20) we have for RHIC (A = 197): a ~ 1.11 

The initial volume Vb in formula (|19fl can be expressed 
across experimental meaning of average number of net 
nucleons N — N in central region of rapidity. For SPS 
the number of net protons in central region of rapidity is 
estimated as 28-29 0. We use the value N - N = 57. 
The number of net protons for RHIC in central region of 
rapidity (at y/s = 130 GeV) is estimated as 8-10 (23, i.e. 
the number of net nucleons is ~ 16-20. We use the value 
N — N ~ 16.3 (for some greater agreement with spectra). 

We have: N-N = V ni(G ,m q , fx q , T, T c ) (from ifT^ . 
Hence the value Tq is: 



[dN 



Hi 



\dy N-N 



1/3 



2tt 2 
12no2 



1/3 



(22) 



We single out here the factor Tq in formula l|13(l for ri\ 
and take into account the values 4^- — 803 and N — N ~ 

dy 

57 for SPS and ^ ~ 1374 andN - N ~ 16.3 for RHIC. 

ay 

In result we have the equation for definition of initial 
temperature Tq: 



Do = a a 



d 2 {G ,m q ,fi q ,T Q ,T c ) 
n 02 {G ,m q , Lt q ,T ,T c ,m s ) 



1/3 



= 1, (23) 



where 



d2(Go,m q ,iJ, q ,To,Tc) = 
= (Fi(Go,m q ,fjL q) To,T c ) — Fi(Go,m a ,- 



■H q >Ta,T c )). 

For calculations of values no 2 we use formulas (I14I17|) 

For SPS we have with accounting of constant factors: 
ao = 1.676 at a = 1.1 and for RHIC ao = 3.045 at a 
= 1.11. We have from here the estimations of initial 
temperature T ~ 175 MeV for SPS and T ~ 219.6 
MeV for RHIC. 

We use here the meaning for temperature of phase 
transition T c = 170 MeV, as indicated by result from 
lattice gauge theory |2^|, and close to T c the tempera- 
ture of chemical freeze-out. At SPS energies we have 
for value fiB — 247 MeV the good approximation for 
different antiparticles and particles ratios. For RHIC 
(y/s = 130 GeV) we use also T c = 170 MeV and fi B = 50 
MeV. That give the ratio p/p w 0.56 close to found be- 
fore experimental value. Later was given |2-l| the ratio 
p/p = 0.6 ± 0.04 ± 0.06. In the work [J were given more 
precise values /is = 46 MeV, T c = 174 MeV and par- 
ticle ratios(at y/s = 130 GeV). However there are for 
the present noticeable statistic and systematic mistake 
in data. One should that not great distinction have weak 
influence for example on spectra baryons and mesons. 

We use also the conservation of entropy and of num- 
ber net nucleons in initial phase at T = To, V = Vb and 
in mixed phase at T = T C ,V = V c (the entropy den- 
sity is s c ): s (T )Vo = s c (T c )V Cl m(T )Vb = ni(T c )V c 



N-N, i.e. V 



The entropy density s c (T c ) correspond to sum (J 
where Tq — > T c . Hence we have relations: 



sq(Tq) _ s c (T c ) 
ni(T ) ni(T c )' 



(24) 



These relations depend also from constant Go- The 
relations l|23|) and (|24l) must be fulfilled simultaneously. 
For SPS (n B =247 MeV, T c = 170 MeV) these rela- 
tions indeed are fulfilled. It can be seen from mentioned 
below values. For example for a = 1.1 at Go — 1.9, 



T 



175 MeV we have : 



so (To) 



Sc(T e ) 



61.4. It 



ni(T ) m(T c ; 

should be noted, that relations l|24(l are fulfilled only for 
coincides values [Xb at initial and mixed phases. Thus, 
in this model the initial temperature To at SPS en- 
ergy exceed weakly the temperature T c of phase tran- 
sition. We will express all physical values in units 
m v = 139 MeV ~ ^ fm" 1 . 
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From formulas (|li p -<|13 p - we find the values of en- 
tropy density for SPS at T and T c : s (T ) ~ 19.43m 3 , 
s c (T c ) ~ 15.31m 3 and also net nucleons density tii(Tq) = 
0.314m 3 , m(T c ) = 0.25m 3 . 

For the number of net nucleons 
(N — N) — ni(T )Vo ~ 57 we find the volume of initial 
plasma phase: Vq ~ 181. 3m~ 3 ~ 518 fm 3 (i.e. tq ~ 3.28 
fm). The initial entropy is: Sq — soVq ~ 3508. We 
have also the volume V c of plasma part of mixed phase, 
taking into account the conservation of entropy and of 
net nucleons: V c = = -^#r ~ 228m" 3 ~ 650 fm 3 

s c (T c ) ni(T c ) 

(i.e. t c ~ 4.1 fm). 

Therefore with accounting of valent u, d quarks we 
have the large initial tq > 1 fm. The suitable calcu- 
lations by formula (|19|1 for To = 1 fm (for SPS) shows 
that for various values of Go the number 57 of net nu- 
cleons correspond fiB — 330 MeV, that contradict to 
experimental data for particles ratios. For correct val- 
ues of /is ~ 200 - 250 MeV the number N - N is too 
little(~ (35 — 45)). The analogous situation there is and 
for RHIC. 

At decrease of interaction the thermal masses of quasi- 
particles m q (T) and m s (T) IjlOjl are diminished at T — > 
T c from above, but they do not correspond to constituent 
masses of quarks in mixed phase. For example, in the pa- 
pers [25l l26j| was shown, that by means of help of masses 
constituent quarks m u = = 363 MeV, m s = 538 
MeV can be found the masses of baryons in octet. For 
SU(3) group the spontaneous breaking of chiral sym- 
metry and appearance of masses leads to appearance of 
octet pseudoscalar light Goldstone (or rather pseudogold- 
stone) states n, k, r\ [23. However if in place of relation 
n^r) * n Q24|) to use relations for massive quarks with ac- 
counting of pseudogoldstone states, then these relations 
do not fulfilled. But these relations will be satisfied, 
if to use decrease of the number of effective degrees of 
freedom of massive quarks. We calculate the entropy 
and net nucleon density s c and n c of massive quarks 
by formulas for ideal gas. For example for m u = 363, 
m a = 538 MeV at [i B = 247 MeV, T c = 170 MeV we 
have s c (T c ) = 14.77m 3 , and n c (T c ) = 0.33m 3 . 

For entropy density s ps of pseudogoldstone states s ps = 
s-k + Sk + s n we find by same way: s ps ~ 3.76m 3 (at T c 
1 70 MeV, m„ ~ 135 MeV). The factor ft of decrease of 
number degrees of freedom can be defined from relation 
g£ ^gf a ^ 61.4. That gives ft ~ 0.685. 

Therefore, it is possible to interpret, that there is ap- 
pearance of massive constituent quarks with effective 
number of degrees of freedom in the presence of octet 
of pseudogoldstone states, while the multigluon state is 
absorbed into constituent quarks masses. 

Let us list appropriate calculations for RHIC 
(y/s = 130 GeV). For example, at a = 1.11 
G = 1.9, T c = 170 MeV we find from ®: 
T ~ 219.6 MeV, that is from here: 

s (T ) 56.44, s c (T c ) 14.8, 

m(T ) = 0.145, m(T c ) = 0.048 in units of m 3 . 



The initial volume V ^ ^ 112.5m" 3 ~ 322.5 fm 3 , 

n i\ 1 o) n 

that is r ~ 2.2 fm (for Tr{R A u) 2 ^ 148 fm 2 ). The initial 
entropy So = sq(Tq)Vo ~ 6350. 

Hence we have for RHIC ^£P4 ~ 389, i#i ~ 308, 

ni{To) ' ni(T c ) ' 

that is the relations (|24J) do not fulfilled. However one 
can assume that at sufficiently high initial temperature 
To > T c it is possible already the appearance (in non- 
perturbative state at T — T c ) of pseudogoldstone with 
entropy density s ps . It gives in fact ^^rf 1 — 387. 
That is we have on this stage of mixed phase for RHIC 
s Vo ~ (s{T c ) + s ps )V c , where V c = ~ 339m- 3 , i.e. 

r c ~ 6.5 fm. 

At high energy (RHIC) the multigluon state is 
also absorbed into massive constituent quarks with 
lesser number of degrees of freedom. Can be found 
the factor ft of decrease of effective number of 
degrees of freedom. We calculate here by anal- 
ogy with SPS the entropy density s c i and net nu- 
cleon density n c \ for massive quarks m u ^ — 363 MeV, 
m s = 538 MeV at /ig = 50 MeV, T c = 170 MeV. We 
find in result: s c i = 14.09m 3 and n cl — 0.0643m 3 
We must have ratio Sclft t Sps ~ 390. We find from 

"clP3 

here ft ~ 0.343. We have effective volume V c i on the 
stage of constituent quarks V c — - ~ 739m_T 3 , or 

s cl P3 T*ps 

V c = N /^r N sf, — 739m~ 3 - is the same. 

Let us consider now the hadronic part of mixed phase. 
We have seen, that entropy density of massive constituent 
quarks is defined by effective number of degrees of free- 
dom and by contribution of octet pseudogoldstone . It 
can be shown that with such effective number of degrees 
of freedom (ft and ft) in hadron part of mixed phase 
appears nucleons and hadrons. In order to find strange 
chemical potential fi s for hadron part of mixed phase, we 
use condition of disappearance of strangeness: 

]T<-5>| = (25) 

i i 

that is: Y\ 9i{F u - F u ) - £j gj (F 2j - F 2j )s = 0. Here 

OO 

f x 2 
Fli = / dx _ , 

gV T 71 T — 1 

Fu = F u (p s -» -fJ, s ), 



F 2j = / dx _ , 

e V x +jv t +1 

F 2 j = F 23 (HB - Sfi a -> -fl B + SfJt, s ). 

For strange baryons we have s = 1 except S(s = 2) 
and Q(s = 3). At temperature T c one should take into 
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account the considerable number of resonances. We took 
into account the strange mesons up to m k = 1820 MeV 
and the strange baryons up to tub — 1940 MeV. In result 
of calculatios we have from fi s ~ 61.44 MeV at 

fi B = 247 MeV and fi a ~ 11.23 MeV at fi B = 50 MeV 
(for T c = 170 MeV). 

At SPS we find the entropy density according to 
calculation by formulas of ideal gas: ~ 5.44m 3 
for nonstrange mesons and ~ 3.26m 3 for strange 
ones (at T c ). For summary entropy density of 
strange baryons and the whole of nonstrange ones 



-p + A + K 



440 



wc find: ~ 5.82m 3 i.e. 



the 



whole entropy is: Sh — 14.52m 3 . For net nucleons 
density with accounting of nucleonic resonances N* the 
calculation gives: nu — 0.402m 3 . Hence the ratio 
is: 2H. ~ 36.2, i.e. ^- < ^. Taking into account 



~ 36.2, i.e. ta- 
riff ' riff 

only nucleons (without of resonances), we have — = 
87.3 > ^. For N + A we have ~ 46.5 



0.1316 

for N + A 



> ^ 

n 



N*. 



520 



v. 



riff 



1535 



the calculation gives 
s H = 13.19m 3 , n H = 0.303m 3 , i.e. ^ ~ 44 and so 
on. Thus always — ^ — , i.e. the relations (12411 do not 
fulfilled. 

However for effective number of nuclcon resonances 
N + A + N^ 440 1520 1535 with lesser number of degrees 
of freedom by factor (3% — 0.685 in the presence of 
octet pseudogoldstones we have the ratio: 13 q^^' 76 = 

eff 

~ 61.4, i.e. the relations (|2"4*j) are fulfilled. We have 
also the volume Vjy = ^jf — 274m~ 3 , or the same 
Vjj ff = =- 274m~ 3 , where V H ff is effective vol- 

riff 

ume of hadron part of mixed phase. We give now the re- 
sult of calculations at RHIC energy. The calculation gives 
for nucleons only (without accounting of resonances): 
530, for n + A this ratio is: ~ 282, for 



10.335 



riff 0.0195 

the whole of nucleon resonance we have for this ratio: 
~ 205 and so on, i.e. — ^ — . But for the whole res- 

' riff ' n c 

onances the calculation gives: sh — 11.85, uh — 0.0585 
in units m 3 . With lesser number of degrees of free- 
dom by value 0$ =0.343(like for massive quarks) we have 

390. Now the effective volume is: 



11.85/33+3.76 



eff 



H 



0.0585ft( — 

812m: 



16.3 

e/7 



or the same ~ 811m„ 



Hence the relations of type l|24|l are fulfilled also for ra- 
tio — in hadron part of mixed but with effective number 

riff r 

of degrees of freedom the same as for constituent quarks. 



IV. THE BARYONS AND MESONS SPECTRA 

We investigate the spectra for SPS and RHIC at 
thermal freeze-out. Apparently the favor results for 
heavy nuclei (for example Pb + Pb) correspond to ther- 
mal freeze-out temperature Tf » 120 MeV |3lll32T| We 
consider at first SPS energy. We find the volume Vf at 
Tf = 120 MeV from conservation of number net nucle- 
ons. The value fi B at thermal freezing we find from 



condition: |(T C ) = §(!», where |(T C ) = e ~W^ 

From here we have \J B ~ 174.35 MeV. From formula 
(|13|l we have net nucleons density at Tf = 120 MeV, 



fjt f B = 174.35 MeV: nf~ N = 7.186 x 10~ 3 m 3 . For 
N net ~ 57 we have V f = N net /n f ~ 7950m" 3 . The 
value n{ at 120 MeV may be found by no condition 
(as the chemical equilibrium already is absent), but from 
ratio: ^{T c ) — ^(Tf) From here we have found with 



accounting of various weak decays: ~ 40 — 42 MeV. 
From relations ^-(T c ) = -^-{Tf) (with accounting 

of decays) we find for SPS the near value: /j{=41 
MeV 331. One can now find by formulas of ideal gas the 
entropy density (at Tf) of nucleons and strange baryons: 

s B ~ 0.125m 3 and for nonstrange and strange mesons up 
to m 4 =1820 MeV: sj l+k ~ 1.3m 3 . Thus the entropy of 
baryons is S B = VfS B ~ 992, and the entropy of mesons 
S[ +k = So — S B ~ 2516. The volume of mesons at ther- 



mal freeze-out is V, 



f 



h+k 



^h + k 

IT — 



1940m 3 We can esti- 



mate the number of ttq: N Vo — V l ^ +k n 7To .For of the value 
n f ~ 0.0807 (calculated without ac- 

counting of weak decays) we have 
N£„ ~ 157. The experiment per unit of rapidity gives 



TO 

we 
N, 



N Vn = 165 ±20 |3J. 
find also 



the 



rf 



0.0149V, 



/ 



By analogy 
number of k + : 
29. Thus in this 



fe + — nk + V h+k — u.u±-±av h+k 
model we have different volume of freeze-out for baryons 
and mesons. The spectra of protons and strange baryons 
one should to calculate with Vf, and of mesons — with 
volume V[ +k . 

The relations of type (|24|l must be fulfilled also 
for thermal freeze-out. Here we have So = Sf = 

s h+k^h+k + s B^f' an< ^ ^ 0r ne ^ nucleons: n{Vo = 
rifVf ~ 57, i.e. we have instead of l|24|) the relation: 
soijh) = 4+ k vLk + 4 ^ 44 + Yl 4 = 61 4 Th th 

ni(To) n f V f n f 

lations of type are fulfilled (but now without of de- 
crease of number of degrees of freedom) And backwards 
from conservation of entropy Sq — Sf one can to find the 
final volume Vf for baryons. 

We calculate the spectra p — p, A, H, SI direct and with 
resonance decays. We take into account the transverse 
flow. We use the hydrodynamic model with linear trans- 
verse velocity profile v±(r/R) = v±(x) — where 
v±_ is the averaged transverse flow velocity. The spectra 
of net protons expressed by formula: 



dNP-P 
m±dm± 



Here cosh(,o) 



Vf 



/Tf 



-m±g p x 



(26) 



sinh(p) cosh(p) 
dxxl (p ± — )Kx (rnj_ -— ) . 



Tr 



sinh(p) 



T f 



vjjx) 



For A, 3, J2 we have factors exp(^ 



- ) , where 
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n = l,2, 3 correspondingly, and for k + - factor 
Vl +k exp(^r-). We take into account also spectra at 

resonance decays, for example: S° — > A7r°(99.54%), 
S_ -» Att- (99.88%), S 1530 ^S7r, S U90 -> Ay, 

Si385 -» Att, n- -> S°tt-. 

At resonance decays we must take into account the 
resonance with transverse flow, and cascade decays also. 

The resonance decays without transverse flow were in- 
vestigated in the paper [3j| . We have calculated here the 
some resonance decays spectra at SPS and RHIC ener- 
gies. For example, spectra of A at decay E° — > A7 is 
calculated thus: 

dN A 

' = (27) 
v=o 



2irm 1 dm I 



T 2 





Here & - branching, g - stat. weight, 
Pi 



+ mi) 2 — m§] [(ms — toi) 2 — m 2 .] 



2m s 

mi = 1115 MeV, 7772 = 0. We have also: 

z (p±) 



f 
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dy 



fu(pj-,V,q,v±)cosh. (y) 



cosh 2 y- f^i- 



(28) 



Here *(p x ) = ln( ^ +K ), ^ = 
The value fn in formula (|2"51) is: 



m s +m 1 — m 2 
2m E 



^2(pj_,a) 



fn(p±,y,q, v±) 



dxx e 



-x cosh (y)^(<?,^± ) 



(29) 



\/NPi^) -x)(x- X!(p±,y)) 
Here r(g,«x) = Vy 1 - (|?«x) 2 , 

/ \ ms(I3imi cosh (j/)=Fp_i_ vABf +p 2 L — cosh 2 (i/) 

a;i,2(p±,yj m i ± cosh i { y } _ p i ± - — • 

At cascade decays, for example E° — > A7, A — > pn, we 
calculate at first the spectra of A, however no at angle 
of 90 degrees (i.e. y = 0), but for the whole y, and then 
this spectra we use for calculation of spectra 2 -JpTdp ± ■ 
In spectra 27rp ^ A rfpx the value exp {-x cosh (y) r(q,v±)} 
we replace now by the value 2K\[x r(q,v±)). Then the 
spectra m d ^ m = f(mx) we substitute to formula: 



m±dm± 

dNP 
2irp±dp±/ y=0 

*(P±) 

dy 



m A bg A 2m A 
Airp* 7rm ± 

dyi 



* /cosh (y A ) 



x (30) 
VifimAy/yi - 1) 



":ei(p_l>I/a) 



V( X 2 - (2/1 - Xx) 



We apply now to RHIC energy (• v /s=130 GeV). At mea- 
surement of antiprotons spectra in collaboration STAR 
were measured the spectra of direct p and also spec- 
tra of p from weak decays (with cascade). The our 
calculation shows that contribution of weak decays is 
~ 30%. We give here the meanings of values, which 
were found for RHIC at T f = 120 MeV (by analogy 



with SPS): t J B =35.3 MeV, /x{=7.93 MeV, 



of 



0.073ml s f h+k 



N—N 



1.287711. From 



conservation of number net nucleons we have the volume 
Vf at Tf = 120 MeV without accounting of weak decays: 
Vf = ~ — 15250m~ 3 . The entropy of baryons in this 



volume is S 



f 



is 5*0 -~ S 



j 



V f s S B 



1.28m 3 



1114. The entropy of mesons 
5236. The mesons entropy density is 
4090m- 3 . 



From here V^ +k 



I 



B ur' 













V ± - 0.44c 








r ** 
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r 
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Figure 1: The baryons spectra in Pb+Pb collisions at SPS 
energy in central region of rapidity (quasiparticle model). 
The A - are direct and with resonance decays: E 1190 — > A7, 

E 1385 _^ ^ H 1315 _^ A?] . -1315 gpectra _ are direct + 

H 1530 — » En + CI" — > H°7r~. The spectra p — p and Qi 672 
- are thermal direct. The averaged transverse flow velocity 
v± — 0.44c. The experimental data are from I29L 1.3911. 

We have found the particles density at Tf. 



0.067m 3 , 



-k+ 



0.0115m 3 . 



This 



gives the estimation of number of particles: 
~ 274, N k + ~ 47. The value V^ +k we use for 
calculation, for example, of spectra k~ and ir + for 
energy RHIC. One can also to show that relation of type 
(|24|l for thermal freeze-out are fulfilled and for RHIC 
energy. 

One can calculate the contribution of weak de- 
cays into net nucleons density at Tf = 120 MeV 



(for RHIC) 33]: 



2.96 x 10" 4 7 



With 
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2.5 3.0 3 J PitQeV/b) 



Figure 2a: The tt° spectra (SPS, Pb+Pb) are direct + decays 
p — > 2tt + f 2 — > 27T (Sx = 0.44c). The dashed line - the direct 
7T°, the dash-dotted line — tv° at decays. We show also the 
summary spectra at v± — 0. The data — from |34l0. 



V? 



Iff 1 



10" 



r 










V -0.44c 








>,., i 


.... i.i.i 


i 



0.0 



02 



0.4 



0.8 1.0 12 ra ± -ro(GsY) 



Figure 2b: The k + spectra (SPS,Pb+Pb) - direct + decay 
K$ so — > kn (v± = 0.44c). We show also the same spectra at 
v± = 0. The data — from [30^. 



account of nucleons from decays we 



lent volume V, 



V f = 16.3 



-N y f 



find equiva- 
From here 



V f ~ 21150m- 3 . 

The spectra of direct and of daughter p at weak de- 
cays we calculate by formulas <|26[l - (13011 with factors 

• • . The contribution 



type Vf exp^s- 



V f exp • 



of weak decays with cascade decays to p spectra gives 
-30%. 

The theoretical p spectra and experimental spectra for 
central bin — 6% (STAR Coll.) at mid rapidity are 
shown in Fig. [3a] In Fig. 13 bl Fig. [3c| we show also the 
direct ir + spectra and with decay p — » irir and the k~ 



(OoV*) 




spectra direct + k. 



x<)2 



kir and experimental data. In 



ni a -m_ (GeV) 



Figure 3a: The spectra p (centrality — 6%) at RHIC energy 
with accounting of weak decays(uj_ = 0.55c). We show also 
the same results with v± = 0.44c and v± — (quasiparticle 
model). The data from f^l. 



Fig.Q]we show the baryons spectra: p—p, A, S, and ft for 
SPS (at Tf) in central interval of rapidity. The tt° and k + 
spectra are shown in Fig.[2aJ Fig. 12 bl together with data. 
One can to see that spectra for RHIC have greater inverse 
slope than for SPS. With increasing mass of particle the 
inverse slope (i.e. the value v±) increases. Of course, one 
should note that obtained spectra with sufficiently large 
value of v± are correct only for not great values of p± . 

Thus, in the present quasiparticle model with isen- 
tropic evolution the normalization of baryons and mesons 
spectra do not defined by free parameteres, but defined 
by conditions in initial state (To, Vq), i.e. by quark-gluon 
plasma. 



V. THE CALCULATION BY THE USE OF 
PERTURB ATIVE THEORY 

We had investigated the physical characteristic of ini- 
tial and mixed phases and some baryons and mesons 
spectra for SPS and RHIC in effective quasiparticle 
model with isentropic evolution. In this section we 



10 




Figure 3b: The ir + spectra (RHIC) - direct + decays p —* tttt 
(v± — 0.48c). It is shown the same spectra with v± =0). 



Figure 3c: The k spectra (RHIC) - direct + decays Kg 92 
kn for v± = 0.52c. The data — from |2ll40| . 



shall consider at first the analogous problem of nu- 
clear collisions by the use of ordinary perturbative 
decomposition of thermodynamic values in powers of 
running coupling a s (T) in the form [41|: a s (T) = 
6n/(UN c - 2N f )\n((2m s + 3T)/A) up to order 0(a s ) 
with QCD parameter A ~ 0.2 GeV. Here the value 
Q 2 in coupling a s (Q 2 ) is considering as the aver- 
age of Q 2 in S-canal at temperature T. For exam- 
ple, the entropy density for gluons in plasma is : 
s g = 32/45tt 2 T 3 (1 - and so on. We calculate here 

also the values T , Vq, So- 

The initial temperature Tq can be find from formula 
analogous J23J), where we have now: 



D (m q , fi q ,T ,m s , X) = 1. 



(31) 



Here the values d,2{m q , a q , To, A) and 
^02(w g , fi q , Tq, m s , A) are extended up to 
order 0(a s ). We find from here for 

a = 1.1, u b = 247 MeV, m q = 0, A = 180 MeV the such 
meanings(for SPS): T = 178 MeV, a s (T ) = 0.454, 
s (T ) ~ 26.4m 3 , s c (T c ) ~ 22.4, m(T ) = 0.475, 
m{T c ) = 0.426 in units m 3 , V = N /nx(T ) ~ 120m" 3 , 
So = soVb = 3180 (where N = (N - iV) ~ 57). 

In this model the entropy of nucleons and strange 
baryons is the same as in quasiparticle model: S 



f 



992, the entropy of mesons S 



f 



f 

h+k 



Sq — S 



f 



1700m" 3 . We have here N 



f 



136, N k 



2200 and 
25. 



h+k — J- • " c ^V — iv *: + 

For RHIC energy we find for a = 1.11, ub = 50 MeV, 
iq = 0, A = 180 MeV the meanings: T ~ 219 MeV, 



m(T c ) = 0.084m 3 , V = iV /ni ~ 111.7, V c ~ 193m" 3 , 
= SoVb ~ 5600. We have here r ~ 2.16 fm, 
t c ~ 3.73 fm. One can show, that in constituent quarks 
phase the relations l|24ll here also are fulfilled only with 
decrease of number of degrees of freedom in the presence 
of pseudogoldstone states. 

We have here by analogy with quasiparticle model 
Vf ~ 15250m 3 , the entropy of baryons S B ~ 1114, 



the entropy of mesons is S 



f 

h+k 



So — S 



f 



V, 



f _ 



h+k 



(s 



ft* = 3500m" 3 

h + k 

the estimation: N n + ~ 243, A^, 



/ 



= 4486 and 
That gives 



a s (T ) = 0.418, s (T ) = 50.11m 3 , m(T ) 



0.145, 



1.28m 3 ). 
fe + - 40. 

In this simple model we have not great decrease 
of normalization of mesons spectra at SPS and RHIC 
(by 12-13 %) in comparison with quasiparticle model, 
i.e. apparently it takes place within the limits of pre- 
cision of measurement. Thus we do not have here no- 
ticeable difference for spectra of particles in comparison 
with quasiparticle model. However this model disagrees 
with SU (3) lattice data in region of phase transition. 
The values s/T 3 , e/T 4 , p/T 4 decrease very weakly in this 
case (For example, at 3T C — > T c the value s/T 3 decrease 
only by value ~ 0.78) The coupling strength a s (T) in- 
crease (although not great) at T — > T c from above (unlike 
effective quasiparticle model). Thus we see, that spectra 
of particles weakly depend on character of phase transi- 
tion QGP — * hadrons. 

But one can show that in this model we have 

too great energy loss AE g of gluon jet at RHIC energy: 
~ 80 — 90%. That is caused in the main by value of 
effective coupling strength. 

Let us consider now the quasiparticle model with phe- 
nomenological parametrization of running coupling G(T) 
© We use here the formulas of type JT0J|-|JT3J| and of 
type H14|) - (|18f) f° r calculation of thermodynamic values 
m q , so and no, but with coupling G(T) (JSJ and with 
C(T,T C )=1. With help of corresponding parameters T s , 
A it is possible well to fit the new lattice data even near of 
T c . For calculation of initial conditions we use the equa- 
tion of type H23f) . For example, for possible parametriza- 
tion 0: A=5.3, T s /T c =0.73 (at N f =S) we find for RHIC 
energy (y / s=130 GeV) the following meanings: To ~ 
250 MeV,G(T ) ~ 2.533 (i.e. a s (T ) ~ 0.509) and also 
a s (T c ) ~ 1.948. We find corresponding values for entropy 
and net nucleons density: so(To) ~ 84.36m 3 , s c (T c ) ~ 
12.96m 3 , m(T ) ~ 0.210m 3 , m(T c ) = 0.0442m 3 . From 
here we have the estimation for initial state: To — (N — 
N)/{ ni {To)iTR 2 Au ) ~ 1.5 fm., and by analogy r c ~ 7.14 
fm. We have for complete entropy: SqVo = 6548, i.e. 
only on ~ 3% the difference from its value in effective 
quasiparticle model. Thus we have practically identical 
baryons and mesons spectra in these two models. 

However we have calculated the jet quenching in 
plasma in this model by analogy with effective quasi- 
particle model 01 • The result is dramatic: the energy 
loss A.E g of gluon jet (for two parametrization of A and 
T s /T c ) four - five times as much than energy E g of itself 
jet !. This contradiction is caused by too great value of 
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running coupling G(T) (JSJ) in this model, which besides 
increase at T — > T c from above. 



VI. CONCLUSION 

In this paper we investigate the isentropic evolution 
of expanding quark-gluon plasma with phase transition 
at first into plasma part and then hadron part of mixed 
phase at temperature T c , and further isentropic expan- 
sion up thermal freeze-out at temperature Tf. For inves- 
tigation of initial plasma phase we use the quasiparticle 
model, which allow apparently to sum up partially the 
perturbative expansion of thermodynamic function of in- 
teracting plasma in powers of coupling constant. In the 
papers jj, |6j was used the phenomenological parameter- 
ization of coupling constant in quasiparticles model in 
accordance with lattice data. However close to T c the 
coupling constant is sufficiently great, therefore in the 
vicinity of the phase transition the perturbative methods 
are not expected to be reliability . 

In the paper was shown, that at such parameter- 
ization the thermal mass of gluon grows in the vicinity 
of the phase transition when T — > T c from above. This 
do not agrees with SU(3) lattice data for Debye mass. 
In the paper ^(j was considered the phenomenological 
model of confinement, where the decrease of thermody- 
namic values at T — + T c from above is caused no increase 
of mass, but by decrease of the number of active degrees 
of freedom (i.e. by modification g g — ► C(T)g g ). It is 
supposed, that behaviour of thermal gluon mass m g (T) 
is similar with behaviour of Debye mass, i.e. it decreases 
with T — > T c from above. This gives the good description 
of SU(3) lattice data. The similar modification is used 
also for quasiparticle model with quarks. 

It is assumed that Debye gluon mass have form 
m D (T) = G(T)T and it has a small gap at T = T c in ac- 
cordance with SU(3) lattice data. It is assumed also that 
the proportionality of type J5J) between mjj and ther- 
mal gluon mass m g remains true in the vicinity of phase 
transition. 

We use such model for description of isentropic evo- 
lution in common case of quasiparticles with /ig ^ 
for calculation of physical characteristic of initial plasma 
phase: To and Vq. In result we have obtained the formula 
(|23[l for definition of these values, where we use the ap- 
proximation (|20[1 for an d experimental mean- 
ing of average number of net nucleons (N — N) in central 
region of rapidity at SPS and RHIC energies. From con- 
servation of entropy and of number of net nucleons follow 
also the relations 1)24(1 . which must be fulfilled simulta- 
neously with J23)- The relations (J23J), depend on 
effective coupling constant Go, which is determined by 
asymptotic value of thermal gluon mass at T = 3T C [l3| . 
For that we use the asymptotic value of running coupling 
ce s {T) @- In result were obtained the initial meaning To, 
V Q at SPS and RHIC energy: for SPS we have T ~ 175 



MeV, V a 518 fm 3 , (i.e. t ~ 3.28 fm) and V c ~ 650 
fm 3 (i.e.T c ~ 4.1 fm), and for RHIC T ~ 219.6 MeV, 
Vo ^ 322 fm 3 and V c ~ 965 fm 3 , i.e. r ~ 2.2 fm, 
t c ~ 6.5 fm. We have also Sq ~ 3508 and ~ 6350 for 
SPS and RHIC correspondingly. 

When the temperature is lowered (T — * T c ), 
more and more gluons form heavy glueballs, and 
the weaker becomes the interaction. This state 
can be interpreted as superposition of nonpertur- 
bative multigluon and weak glueball exchange [To| . 
The interaction between the quarks also becomes 
weaker. It is possible to interpret, that there is 
appearance of massive constituent quarks, while the 
multigluon state is absorbed into constituent quarks 
masses. 

Further we show, that massive constituent quarks ap- 
pears with decrease of number of degrees of freedom in 
the presence of octet of pseudogoldstone states. That 
follow from conservation of entropy and of number of 
net nucleons. The hadron part of mixed phase appears 
also with decrease of number of degrees of freedom - 
the same as for constituent quarks. However at ther- 
mal freeze-out (at T = Tf) hadrons appears already with 
normal number of degrees of freedom. 

Thus we have at SPS the short plasma phase. With 
above-mentioned parameters we have calculated in ef- 
fective quasiparticle model the some typical baryons and 
mesons spectra with account of resonance decays. We use 
for heavy nuclei the temperature of thermal freeze-out 
T f = 120 MeV (as in "standard" model |MI13 The nor- 
malization of baryons and mesons spectra in the present 
model is no free parameter, but defined by initial condi- 
tion in quark-gluon plasma. 

One should note that there are the models for descrip- 
tion of particle production with simultaneous chemical 
and thermal freeze-out 36], where 2-3 free parameters 
describe well the spectra at RHIC energy. In the present 
work we use "standard" scenario with two freeze-out - 
chemical and thermal. 

The main object of present work is no study the whole 
spectra at SPS and RHIC, but the study of dependence of 
spectra from initial state and from phase transition stage. 
For this aim we consider in Sec0the analogous isentropic 
problem by the use of ordinary perturbative theory up to 
order 0(a s (T)), where a s (T) is running coupling. We do 
not find in this approximation of noticeable difference for 
particles spectra in comparison with effective quasipar- 
ticle model, although this model disagrees with SU(3) 
lattice data in region of phase transition (unlike quasi- 
particle model). However it can be shown that in 
this model we have the too great energy loss of gluon 
and quark jets in plasma (~ 80 — 90%), that contradict 
to suppression of 7r°-spectra, reported by PHENIX (Ijj . 

In this section we consider also the quasiparticle 
model with phenomenological parametrization of cou- 
pling G(T), which gives good fit of new lattice data. We 
investigate here also the initial condition in plasma at 
RHIC energy (on the basic of equation of tvpe (|23|l . We 
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find that baryons and mesons spectra in this model prac- 
tically do not differ from spectra in effective quasiparticle 
model. 

However the situation with energy loss in plasma here 
yet more dramatic - the energy loss of gluon jet exceed 
the energy of itself jet. Thus apparently we have cor- 
rect quantitative description of energy loss of gluon and 
quark jets in plasma only in considering here effective 
quasiparticle model. 

But of course the question about nonperturbative so- 
lution for thermodynamic values in quark-gluon plasma 
remains open (for example, for thermal mass of quasipar- 
ticles and running coupling and so on). Unfortunately 
in the vicinity of the phase transition is used the phe- 
nomcnological parameterization for these values. 

In this model we do not consider also the intrinsic 
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